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It has long been believed that equilibrium random-field 
Ising model (RFIM) critical scattering studies are not feasi- 
ble in dilute antiferromagnets close to and below T C (H) be- 
cause of severe non-equilibrium effects. The high magnetic 
concentration Ising antiferromagnet Feo.93Zno.07F2, however, 
does provide equilibrium behavior. We have employed scaling 
techniques to extract the universal equilibrium scattering line 
shape, critical exponents v = 0.87 ±0.07 and r\ — 0.20 ±0.05, 
and amplitude ratios of this RFIM system. 



Models for the statistical physics of phase transitions 
can be experimentally tested by measuring the universal 
critical parameters and comparing them with theoreti- 
cal predictions. This is crucial for the understanding of 
difficult problems such as the random-field Ising model 
(RFIM). Although it is well known that a phase transi- 
tion occurs J!]] for the d = 3 RFIM, the characterization 
of the equilibrium critical behavior in its experimental 
realization, anisotropic randomly dilute antiferromagnets 
such as Fe x Zn\— x F% in applied uniform fields has not 
been possible near to and below the transition, T C (H), 
despite years of intense experimental investigation. This 
is primarily a result of the severe hysteresis in the scat- 
tering line shapes near to and below T C (H). Although 
long-range order is observed after cooling in zero field 
(ZFC) and heating through the transition, upon cool- 
ing through the transition in the field (FC) long-range 
antiferromagnetic order is never achieved, even at low 
temperatures. In either case, the line shapes are difficult 
to interpret. It was widely thought for many years that 
these nonequilibrium effects were unavoidable. 

Since random-field effects increase with the dilution, as 
well as the strength of the applied magnetic field || , most 
early studies focused on magnetic concentrations x < 0.8 
in order to achieve, for reasonable H, a suitable range 
of reduced temperature, t — (T — T C (H))/T C (H), over 
which asymptotic critical behavior might be observed. 
However, nonequilibrium effects exist below the equilib- 
rium temperature [0 T eq (H) which lies just above T C (H). 
Much effort was directed to understanding the behavior 
of the domains, both near to and below T C (H). Neu- 
tron scattering studies of the Bragg intensity in the bulk 
crystal Fcq, 46^^0.54-^2 and the epitaxial film || of 



FeQ.52Zna.4sF2 showed anomalous behavior even under 
ZFC caused by domain formation. It was shown that, at 
low T, RFIM dynamics are dominated by domains 
The dominance of domain effects near T C (H) was evident 
in ac susceptibility measurements || . From these studies 
over several years along with Monte-Carlo (MC) simula- 
tions [^o|JTl| l it became clear that the severe nonequilib- 
rium effects were primarily due to the large number of 
vacancies that are well connected and thus allow the for- 
mation of domain walls with insignificant energy cost. 

Recently, we demonstrated Q| that the metastability 
problem can, in fact, be overcome by employing a mag- 
netically concentrated crystal, Feo. 93^710.07^2, in which 
the domain walls must break many magnetic bonds to 
form since the vacancies are well separated As a 

consequence, domains do not form under either FC or 
ZFC and the line shapes from both processes are found 
to be identical for all T and, hence, represent equilibrium 
behavior. Of course, the small dilution necessitates the 
use of large fields to probe the RFIM critical regime, so 
we used the highest field available with the present appa- 
ratus, H = 7 T. Attempts to analyze the equilibrium crit- 
ical scattering in a preliminary study [ fl2| fa iled because 
the sample had a concentration gradient KM that limited 
analysis to \t\ > 10~ 3 . Hence, we were only able to at- 
tempt an analysis using mean-field (MF) line shapes that 
yield values for the inverse correlation, n, and staggered 
susceptibility, \ s , that do not, however, follow power-law 
behavior in \t\ below T C (H). We have overcome these 
difficulties by using a thin crystal of thickness 0.44 mm, 
less than 1/10 the original thickness, cut with its plane 
perpendicular to the concentration gradient. This has 
made it possible to extend the scattering measurements 
to \t\ > 10~ 4 , greatly increasing the accuracy of the crit- 
ical behavior analysis and, by making a few simple as- 
sumptions about the scaling properties of the scattering 
line shape as described below, we have been able to go 
beyond the inaccurate MF line shapes traditionally used 
in this type of analysis. With the use of the scaling tech- 
niques and the additional data from the thin crystal, we 
have finally achieved a consistent analysis of the critical 
behavior above and below T C (H) for the RFIM. 

Neutron scattering measurements were made at the 
Oak Ridge National Laboratory High Flux Isotope Re- 
actor using a two-axis spectrometer configuration. We 
used the (0 2) reflection of pyrolytic graphite (PG) at 
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an energy of 14.7 meV to monochromate the beam. We 
mainly employed two different collimation configurations. 
The lower resolution, primarily used for the large sam- 
ple, is with 70 min of arc before the monochromator, 20 
before the sample and 20 after the sample. Primarily for 
the thin sample, we made scans with 10 min of arc before 
and after the sample. PG filters were used to eliminate 
higher-order scattering. The carbon thermometry scale 
was calibrated to agree with recent specific heat results 
p5[ for the H = transition. The field dependence of the 
thermometry was also calibrated. All scans used in this 
report are transverse about the (10 0) antiferromagnetic 
Bragg point. 

The observed line shape is given by S(q) convoluted 
with the instrumental resolution. In MF 
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where n = K ± |t| iy , + and — are for t > and t < 0, 
respectively and M s is the staggered magnetization. The 
first term represents Xs (q) and the second two the discon- 
nected susceptibility, Xs dls {q)- In principle, the ampli- 
tudes A and B are constant in T and the same above and 
below T C (H), although in practice this is not true since 
real systems do not follow MF behavior. For translation- 
ally invariant systems, B — 0. For real systems MF is in- 
adequate since we must have the scaling behavior Xs (q) oc 
kV 2 f(q/ K ) with the limits Xs(q) tx k''~ 2 /(1 + q 2 /n 2 ) for 
\q\ << k and Xs(q) of q v ~ 2 f° r \q\ >> K - For the d = 3 
pure transition, r\ w 0.03. Hence, the MF Lorentzian line 
shape is adequate except very close to T C (H), as explic- 
itly shown EJ| for FeFi where the corrections to MF are 
important for \t\ < 10~ 3 , particularly for T < T C {H). 
For smaller the t < Tarko-Fisher EJ] (TF), 
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approximants, where cf>, </>', 0", and a are fixed, and ip = 
l + l/2r/0 2 and ip' = l + l/2?70 /2 -|-cr(>' 2 -</>" 2 ), were found 
to be excellent for the analysis and yielded exponents in 
superb agreement with theory EJ . 

The dilute case with H = corresponds to the random- 
exchange Ising model for which the value of rj is similarly 
small and the MF line shape proved adequate for \t\ > 
10~ 3 . For x = 0.46, the scattering results EJ| v = 0.69 ± 
0.03 and 7 = 1.33 ± 0.02 agree remarkably well with 
recent MC results [§0| v = 0.684 ± 0.005 and 7 = 1.34 ± 
0.01. 

In a previous scattering work Ej ] , we found two impor- 
tant results using Eq. 1 to fit the data for Fe^.^Zn^mFi 



for H = 7 T. First, the results for k and Xs are indepen- 
dent of whether the FC or ZFC procedures were used, i.e. 
the line shapes are equilibrium ones at all T. Second, the 
values for k and Xs obtained from the line shape anal- 
ysis could be fit to a power law in \t\ for T > T C (H), 
but not T < T C {H). For T > T C (H), v = 0.90 ± 0.01, 
7 = 1.72 ± 0.02 and 7 = 3.0 ± 0.1 were obtained. The 
lack of power law behavior for k and Xs below T C (H) 
indicates that the MF line shapes are inadequate, partic- 
ularly for T < T C (H). This is not surprising, since the 
predicted d — 3 RFIM exponent rj « 0.5 is quite large ffl 
and a large deviation from MF behavior is therefore to 
be expected, particularly below T C (H), at |t| much larger 
than in the pure FeFi case. 

Since we lack a theoretically predicted line shape for 
data analysis, we instead invoke the scaling form EJ 

S(q) = A ± ^- 2 f(q/ K ) + B ± A ±2 ^- 4 g(q/n) (4) 

for \q\ > where 7 = v(2 — rf) and 7 = v(A — fj). This 
expression is still too complicated to use for scaling of the 
data, particularly since the scaling in this unusual case 
involves two perhaps independent functions. To proceed, 
we make two assumptions strongly motivated by results 
of Monte Carlo simulations |2l],^2| , high temperature ex- 
pansions p3[ , and previous experiments |Q at x = 0.6, 
namely fj = 2rj and g{q/n) = f 2 {q/^)- With these as- 
sumptions, we have 



S(q) 



f(q/ K )(l + B ± A ± ^- 2 f(q/K)) (5) 



for q I > 0. Finally, we employ the TF and FB approxi- 
mants for f(q/n) except that we let </>, 0', 0", and a be 
fitting variables. The line shape analysis includes a small 
fixed constant term and one that is linear in q to account 
for background counts. 

For H = 0, the random-exchange Ising model case, we 
set B — 0. Folding the resolution corrections into the 
scaling line shape, we fit all of the data simultaneously 
over the range 1.14 x 10~ 4 < \t\ < 10~ 2 for the small 
sample and 1.15 x 10~ 3 < |t| < 10~ 2 for the large one. 
We obtain the critical parameters v = 0.70±0.02 and 7 = 
1.34 ± 0.06, which are in good agreement with previous 
results using the MF equation cited above. 

The same procedure is followed for H = 7 T, the RFIM 
case, where we now let B vary. The fitted parameters are 
given in Table I for \t\ < 10~ 2 and \t\ < 3 x 10~ 3 . The 
lower limits for \t\ are the same as for H — 0. In Fig. 1 we 
plot, versus q/n, the intensity data after subtracting the 
background, deconvoluting with the instrumental reso- 
lution, and dividing by A ± k'^ 2 (1 + B ± A ± k^ 2 f(q/ '«)). 
The collapse of the data onto f(q/n) is excellent. We 
find no evidence of systematic deviations of the data 
from the scaling function in any of the scans used in 
the analysis. The collapse of the data onto a single scal- 
ing function can only occur if the critical exponents are 
chosen appropriately. By the quality of the scaling, we 
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have demonstrated that by moving beyond the simple 
MF line shape, it is possible to fit the RFIM data for high 
magnetic concentration to obtain the critical parameters 
and a good representation of the actual scattering line 
shape. Note that these values, obtained for data above 
and below T C (H), agree rather well with the exponents 
obtained with the MF expression for this system above 
T C (H). The curves in Fig. 1, given by the TF and FB 
parameters in Table 1, represent the experimental RFIM 
scaling functions that can be used to test future theoret- 
ical and simulation results. 

Numerous simulation and series expansion studies 
yield critical exponents for the RFIM. The recent Monte 
Carlo results ofRieger (fl) are v = 1.1±0.2, 7 = 1.7±0.2, 
n = 0.5±0.05, 7 = 3.3±0.6, fj = 1.03±.05, = 0.0±.005, 
and a = —0.5 ± 0.2. Similarly, Hartmann and Nowak 
p2[, using exact ground state simulation techniques, find 
v = 1.14 ± 0.10, (3 = 0.02 ± 0.01 and 7 = 3.4 ± 0.4. 
Note that the values for v and 7 agree fairly well with 
the results of the neutron scattering experiments (using 
7 = i/(2 — 77)). However, 77 is larger than the value from 
the neutron study. The value of f3 has not been accu- 
rately determined experimentally. The largest discrep- 
ancy p5{ is with a which appears to be close to zero 
since the experiments at all concentrations exhibit sym- 
metric, nearly logarithmic peaks in the specific hea t p| . 
This disagreement has not yet been fully explained |22|. 

The equilibrium exponents obtained in the present ex- 
periments can be compared with previous results at lower 
concentrations. Belanger, King and Jaccarino p3] ob- 
tained v = l.Q ± 0.15, 7 = 1.75 ± 0.20, 7 = 3.5 ± 0.3 
and rj « 1/4 using a sample of Feo. qZuoaF^ and the MF 
equation. Only data above the equilibrium boundary 
T eq (H), were used. These results, are in good agree- 
ment with the ones obtained using the MF line shapes in 
Feo.93Zno.07F2- This agreement indicates that the cor- 
rections to MF are not very significant above T C (H) but 
are very large for T < T C (H). 

In contrast, Feng, et al. |^6| obtain v = 1.5 ± 0.3, 
7 = 2.6 ± 0.5 and 7 = 5.7 ± 1.0 for Fe a , 5 Zn . 5 F 2 . This 
discrepancy, however, may be the result of the analysis of 
data which may lie below T eq {H), where hysteresis pre- 
vents a study of the equilibrium critical behavior. This 
seems apparent Ej]] when comparing these authors' re- 
sults with previous studies ]32|| of Fe^Zn\- x F%. The 
present study at high magnetic concentration does away 
with metastability effects entirely for the critical scat- 
tering line shape analyses. This has the particular ad- 
vantage of allowing fits to the data on both sides of the 
transition and much closer to the transition, greatly nar- 
rowing possible interpretations of the data and yielding 
much more reliable critical behavior parameters. 

In summary, we have characterized the critical scatter- 
ing below T C (H) in the high magnetic concentration crys- 
tal FeomZno.o7Zri2. This has not been possible until 
now because of the domain formation present in the more 



diluted systems previously studied. We have shown that 
the line shapes are adequately described by using sim- 
ple scaling hypotheses. The possibility still exists that 
the exponents obtained in Feo. 97^^0.07^2 are effective 
ones that may change somewhat with fields much larger 
than H = 7 T. Future efforts will be directed to scatter- 
ing measurements to determine v, 7 and (3 and birefrin- 
gence and Faraday rotation measurements to determine 
the specific heat critical behavior for H >> 7 T at high 
magnetic concentration. To avoid extinction effects, the 
Bragg scattering exponent determination requires future 
experiments either using neutron scattering on a film || 
or using magnetic x-ray surface scattering. 
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FIG. 1. Scaled neutron scattering data taken at different 
temperatures at H — 7 T collapsed onto the universal func- 
tion /(g/fc). The scatter in the small sample data is larger 
due to smaller number of counts obtained in the thin sample. 
The fit was made for \t\ < 0.01. 



parameter 


\t\ < 10 -2 


\t\ < 3 x 10 -3 


T c (fixed) 


70.61^ 


70.61A- 


n 


0.20 ± 0.05 


0.16 ±0.06 


V 


0.88 ±0.05 


0.87 ±0.07 


A+ 


10.0 ± 0.2 


9.21 ±0.3 


A~ 


6.15 ±0.14 


4.45 ±0.15 


K + 


1.13 ±0.04 rlu 


0.95 ±0.17 rlu 


K (> 


3.24 ±0.11 rlu 


2.78 ± 0.5 rlu 


B+ 


(4.7 ±0.1) x 10~ 5 


(3.00 ±0.13) x 10 


B~ 


(4.0 ± 0.3) x 10" 5 


(8.0 ± 1.0) x 10" 


a 


0.67 ±0.5 


0.86 ±0.6 


4> 


0.16 ±0.04 


0.08 ±0.01 


4>' 


0.39 ±0.25 


0.36 ±0.3 


0" 


0.31 ±0.25 


0.26 ±0.2 


X 2 


3.07 


2.3 


No. pts. 


2444 


1000 



TABLE I. Values found from fits at H = 7 T. Data for the 
large sample were fit for \t\ > 1.15 x 10 -3 and for the small 
sample were fit for |t| > 1.14 x 10~ 4 . k is given in reciprocal 
lattice units (rlu). 
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